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We consider the p-spin spherical spin-glass model in the presence of an external magnetic field as a 
general example of a mean-field system where a one step replica symmetry breaking (1-RSB) occurs. 
In this context we compute the complexity of the Thouless- Anderson-Palmer states, performing a 
quenched computation. We find what is the general connection between this method and the 
standard static 1-RSB one, formulating a clear mapping between the parameters used in the two 
different calculations. We also perform a dynamical analysis of the model, by which we confirm the 
validity of our results. 
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I. INTRODUCTION 



Two different sets of results can be obtained when studying a statistical model. The first set is given by a static 
(thermodynamical) analysis of the system, while the second one is obtained from a dynamical approach (eg., Langevin 
dynamics). In the case of disordered models as spin-glasses, these two sets of results are not in general connected in 
a straightforward way. The reason for this is that disordered models always exhibit an out of equilibrium dynamical 
behavior, which is believed to be related to the presence of many metastable states. 

However, there is a striking correspondence between certain types of dynamical scenarios and their specific static 
counterparts. In other words, some particular dynamical results are always associated with the same set of static 
results and vice-versa. For instance, in the context of the static replica symmetry breaking (RSB) scheme for mean- 
field spin glasses we always find a dynamical asymptotic energy larger than the static one in models which are 
solved by a one step solution (1-RSB), while the two energies are the same in models statically solved by a full RSB 
solution [^|-^| • This fact suggests that there must be some kind of underlying general physical explanation for these 
connections. 

There is a third method which can be used in the study of disordered systems, which has the advantage to catch 
some of the results both of the statics and of the dynamics. This method is purely entropic, that is it just deals 
with the number of states of the system, disregarding their thermodynamical Boltzmann weights. In this context, 
one simply computes the number of local minima of some mean-field free energy, function of the local magnetizations 
of the system, and obtains in this way the spectrum of the states as a function of their free energy density 
This approach gives information both on the statics and the dynamics of the system since it deals with the whole 
underlying free energy landscape. We believe that in this context it should be possible to find an explanation for the 
relation we stated above between statics and dynamics. 

As a first step in this direction it is necessary to understand how general and feasible this entropic approach is. 
Up to now it has given correct results only in one specific mean-field model with 1-RSB solution, the p-spin spherical 
spin-glass in zero magnetic field |§,||J8). In this model there is a band of free energy densities where the number of 
states i sexponentially large in the size of the system. It has been found that the static and the dynamic free energies 
(the last being defined as the free energy of the states reached by the dynamics) are equal to the lowest and highest 
edges of this band, respectively. We believe this to be true in any 1-RSB model. Yet, in order to find out how general 
this result is, it is not sufficient to consider the p-spin model in zero field. Indeed, due to the absence of the field, this 
model has some particular features, which make it a special case in the 1-RSB class. It is homogeneous and therefore 
presents no chaoticity of the states with respect to the temperature j^] . The static and dynamical transitions are both 
discontinuous, which is not a general property of 1-RSB models. Finally, the dynamical off-equilibrium behavior of 
the p-spin model with zero field turns out to be much simpler than in the general case Q|. However, this model can 
be easily generalized by the introduction of an external magnetic field. It then becomes a very general paradigm of 
1-RSB models, since as a function of temperature and field it displays all the possible characteristics of these systems 

In the present paper we extend the entropic calculation mentioned above to the p-spin spherical model in presence 
of a field. In this context we are able to state in the most general way which are the connections with the standard 
static and dynamical approaches. In particular, we find a one to one mapping between the 1-RSB static parameters 
and the ones of the entropic calculation, giving a physical interpretation of this relation, until now rather unclear. 
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The paper is organized in the following way. In Section |l| we define the model and compute the average logarithm 
of the number of states, known as complexity. In the next Section we discuss the results and find the relation with 
the static 1-RSB calculation. In Section [V we study the behavior of the complexity as a function of temperature and 
magnetic field. In Section ^ we perform a dynamical analysis of the model and compare it with the results coming 
from the complexity. Finally, we state our conclusions in Section VI . 



II. MODEL AND COMPLEXITY 



The p-spin spherical model is defined by the Hamiltonian, 

H{s) = - Jt 1 ...t p s il . . .s lp - h^Si. (1) 

i 1 <...<i p i 

The spins s are real variables satisfying the spherical constraint ^ sf = N, where N is the size of the system. 
The couplings J are Gaussian variables with zero mean and variance p\/2N p ~ 1 and h is an external magnetic field 
|^o| , p] J^1. In the context of the Thouless- Anderson-Palmer (TAP) approach Q, it is possible to formulate a set of 
mean field equations for the local magnetizations rrii = {si). A mean-field free energy density /tap, function of the 
magnetizations rrii, was introduced for this model in jjj), 

/TAP = -^r J h-ip m ii ■■■™> ip ~ Yl mi ~ ^g lo g( 1 - ( l)+9{q), (2) 

ii<...<i p i 

where q = N^ 1 ^ mf is the self-overlap related to the magnetization m and of a) = —4 \{p — l)q p — pq 9 ^ 1 + l] is 
the Onsager reaction term (for a derivation of the TAP free energy see also fl2|,p|). The minimization of (Q) with 
respect to the {rrii} gives the TAP equations of the system, 



T k (m) = -p Y Jk,i 2 ...i p m l2 ...m ip -h + 2m k ( — ^— - + g'(q)j = 0. (3) 



In the low temperature phase these equations admit many possible solutions, corresponding to different stationary 
points of the mean field free energy (||) . The minima among all the stationary points can be identified with stable 
and metastable states of the system. 

The complexity £(/) is defined in the following way, 



£(/) = lim - logAf(/), (4) 

where A/"(/) is the number of local minima of the free energy (|2|), that is, the number of states of the system with 
free energy density /. We average the logarithm of N since we expect this to be the extensive quantity. In order to 
perform this average it is necessary to introduce replicas. However, when the external field h is set equal to zero, the 
correct ansatz for the overlap matrix turns out to be symmetric and diagonal. This is equivalent to averaging directly 
the number M of solutions (annealed average) || . The physical reason for this is that when h = there is no preferred 
direction in the phase space sphere, thus the typical states are orthogonal to each other and their mutual overlap is 
zero. On the other hand, when h ^ there is a migration of the states towards the direction of the magnetic field 
and their mutual overlap is different from zero. In this case the quenched average has to be performed. Anyhow, as 
it is shown below, due to the particular nature of the calculation, the correct ansatz for the overlap matrix is simpler 
than the 1-RSB used in the statics Q . 

In order to compute £ we use the replica trick, 



log N = lim - log TT Af a , (5) 

a=l 

where each M a is given by 

N 

Af a (f)= Vm a Y[6(T k (m a ))\detA(m a )\6(f TAP (m a )-f), (6) 



fc=i 
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and Aki(m) — 0fe<3i/xAp("i) is the Hessian of the TAP free energy evaluated in a particular solution m. In what 
follows the symbol T>x refers to the integration over all the site variables, Dx = dx\ . . . dxpf. Let us introduce the 
following Grassmann representation for the determinant, 



det„4(m a ) = / V^ a Vi) a exp \-iJ> a A{m a ^ a ] 



and a bosonic representation for the delta functions, 

N 



f[S(T k (m a )) = f PA Q exp(-A a T(m a )), 
k=i J 

S (hAp(m a ) - /) = / dcu a exp [~uj a (f TAP (m a ) - /)] 



(7) 



(8) 



(9) 



where the integrals over the variables A and ui are on the imaginary axis. The sums over site indices are always 
understood whenever site dependent quantities are mutually multiplied. In the rest of the calculation we shall 
disregard the modulus of the determinant in (^) . This approximation is safe for zero magnetic field, as long as we are 
counting stationary points with a given fixed free energy density [jl3| . Indeed, it can be shown that in this way only 
the minima of /tap are actually taken into consideration [fill. We assume this to hold also with h 0. We then have, 



£(/) = lirn— ^rlog / Vm V A "Dtp "Dip duj exp [— Sj(m, A, ip, ip, <*>)] . 
where the action Sj is given by 

n 

Sj = [A Q T(m Q ) + $ a A(m a )i/j a + iYa, a (/ T Ap(m a ) - /)] . 



(10) 



(11) 



The average over the disorder generates couplings between variables with different replica indices. It is convenient to 
express these terms by means of the following overlap matrices, which can be introduced in the usual way PJq,|l3|] , 



NS ab = m a m\ NL ab = X a X b , NR ab = m a \\ NT ab = 
In terms of the overlap matrices ( p^ ) we have, 

£(/) = lim — log / VSVLVRVTdu exp \—NS(S, L, R, T, u>)] , 
nN J 



(12) 



(13) 



where the effective action is given by 

• * = z E L ^ x + ^r 11 E - Z 

i 



4 ao | 

E L ^ + E 2 



ab 

if 

2 



a b 

g'(s aa ) + 



(T aa Raa) log det T 



2/3(1 - S a 

\ log det S + \ log det (R T S~ 1 R - L) + ~uj RabS^ 1 + h 2 u £ R ab 



+nfu> — lo 



g(S aa ) - —\og(l - S aa ) 



: E 5afc 



(14) 



We have assumed Lo a = to since it depends only on one replica index. As usual, we compute the integral ( |l3| ) by means 
of a saddle point approximation, so that, 



E(/) = lim -Ext* [S(j, 

n— >0 n 



(15) 



where <f> stands for all the variational parameters, S,L,R,T and lo. 

In order to simplify the saddle point equations for the overlap matrices we make use of the fact that action (|ll]) is 
invariant under the following Becchi-Rouet-Stora-Tyutin (BRST) transformation |l4]-|l(| , 
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m?-mf + e<, #->#- e A?, A«^A«- W " £ C <->V?, w a ^ u/ 1 (16) 

where e is a constant Grassmann parameter. As a consequence, for each operator O, which is a function of the 
variables m", A°, V> a , ip a and w a , we have that (SO) — 0, where SO is the variation of O under ( ffq ) and the brackets 
indicate an average over the measure defined by the action Sj. If we consider the two cases O = m?ip a and O = X b ip a , 
we immediately get the equations (ip a, i(j b ) = —(m b \ a ) and (u)i[) a ip b ) = (X a X b ). From the definitions ( |l2|) we then 
obtain the following relations for the saddle point values of the overlap matrices, 

Rab = T^, L a b = —LL>T a b (17) 

which simplify a lot the calculation. At this point we have to choose an ansatz for the overlap matrices in order to 
solve the saddle point equations. Let us consider first of all the matrix S, which has an explicit physical meaning. 
From the definition jl^ ) we see that the diagonal element S aa corresponds to the self-overlap of an individual TAP 
solution. On the other hand the off-diagonal elements S a b correspond to mutual overlaps between different TAP 
solutions. Given this we assume for the matrix S the simplest form consistent with this interpretation, that is a 
symmetric matrix, 

S a b = (si - s )S ab + s . (18) 

In the saddle point, s± and sq represent respectively the self-overlap and the mutual overlap of the states with free 
energy density /. As discussed above, in absence of the magnetic field s vanishes, the matrix S is diagonal and the 
quenched average coincides with the annealed one. On the other hand in the presence of the field we expect a value 
so ^ 0. Consistently with the ansatz used for S we set, 

T ab = (ti - t )S ab + 1 . (19) 

Using the BRST relations (0) and the above form for the matrices S and T, we can reduce the saddle point equations 
to the following five coupled equations, 



y-t% +hj+ i So -^(ujso + to), 



» . (f » r ' + ^) + ^ - ^ - \ + + w 







(|^ + /l 2) + _^___ (a;S0+ i Q)) 



-{■ --Ins 

~ + ^2^^ s °) + f 9ist) + ± log(l - Sl ), (20) 

where z = [u>(si — sq) +t\ — to] and g(s) is the Onsager term of (||). It is easy to verify that for h — these equations 
give the standard result of § . 



III. CONNECTION WITH THE STATICS 



We are now in the condition to compute the complexity E in every point of the plane (T, h) and in particular in 
the regime where the statics of the model display a 1-RSB solution (low temperatures and fields) As we will show 
below, our results clarify the relation between the present entropic approach and the usual static one. 

First we note that equations (^) give, after some algebra, the relations, 

to=0 , h=(3{l- Sl ). (21) 

The physical interpretation of these equations is very simple. The saddle point value of the matrix T a b is related to 
the expectation value (-0 a ^ b ) (see Eq. (|l2|)) and from (^) it is clear that this expectation value is nothing else than 
the average of the matrix A^ 1 , which is by definition the inverse of the Hessian of the TAP free energy. Therefore, 
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the parameter t\ can be interpreted as the inverse curvature of the free energy in a typical state with free energy /. 
On the other hand, since s\ is the self-overlap of this state, (3(1 — Si) is the the magnetic susceptibility \- I n this way 
equation (|2l]) gives the expected static relation between the fluctuation t\ and the dissipation x of an equilibrium 
system. In other words, the size (1 — s±) of a state has two different contributions, a thermodynamic one proportional 
to the temperature and a geometric one proportional to the inverse curvature of the mean field free energy. Equation 
( pl| ) has not been recognized before in the context of this calculation H|J , simply because for h — the homogeneity 
of the p-spin model permits to compute everything at zero temperature where s% = 1, so that the above relation 
remains hidden. 

Due to equation (^l|) we are left with only three parameters, s±, So and lo. As we have seen, the physical interpre- 
tation of si and sq is straightforward, while this is still not the case for lo. In order to better understand what is the 
role of lo we note that in the saddle point we have (see equations (111) and (ph), 



df 

If we call fo the ground state free energy, defined by E(/o) = 0, the number of states having extensive free energy 
Nf = Nf + 6f, is 

Af(Sf) ~ e E ' (/oW = e " (/oW , df ~ 0(1). (23) 

This equation shows that w(/o) plays the same role as the static parameter x, which represents the breaking point 
in the 1-RSB scheme |r^ , |l8| ,|j| . Given this, it is clear that, for each value of (T, h), the 1-RSB computation and the 
present one must map one into the other with the following rule: 

fo = /rsb => si(/ )=gi, s (f Q ) = qo, Lo(f )=l3x, (24) 

where /rsb is the static 1-RSB equilibrium free energy and q\, qo and x are the 1-RSB parameters of the static 
computation (respectively, self-overlap, overlap and breaking point) It is possible to verify analytically that the 
relations ( p4| ) are fulfilled by our saddle point equations: the complexity is zero at the static free energy /rsb an d 
here the equations for the parameters s\, so and to reduce exactly to the static equations for q\, qo and {3x. 

In this way we have shown that, despite the symmetric form of the overlap matrices used in our calculation (see 
Eqs. (|l8|) and (|l9|)), the complexity describes the statics of the model with the same degree of accuracy as the 1-RSB 
calculation. This point has always been rather mysterious, since it was not evident why the two approaches should 
give the same results. Now the explanation is clear. At the ground state free energy density fo the two calculations 
simply write in different ways the same quantities. In the RSB approach the self-overlap q\ of the states is introduced, 
while for the complexity it is natural to deal with the curvature of the states, because of the presence of the Hessian. 
Yet, these two quantities are, as we have seen, trivially related. Moreover, the breaking point x, which is typical of 
the RSB calculations, has its counterpart in the parameter lo, which has a very simple interpretation in the context 
of the complexity, being just the derivative of E with respect to the free energy. Finally, the overlap matrix in the 
two cases has a slightly different physical meaning and this is why different ansatz are taken for them: in the RSB 
approach this matrix refers to the overlap between configurations, while in the present context it refers to overlap 
between states. Due to this, the present approach has one less step of replica symmetry breaking, compared with 
the standard static one. Clearly, the price we have to pay for having a symmetric overlap matrix is the knowledge of 
the explicit form of the TAP free energy. Interestingly enough, the parameter x drops from the symmetric overlap 
matrix, to reappear in disguise as lo. We believe this scenario to be true in each model exactly solved by a 1-RSB 
static computation. 

We stress that what is stated above only holds for f = fo, where the complexity is zero. On the other hand, 
increasing / provides information on the whole spectrum of the states. 

Finally, we would like to recall a completely different method by which the complexity can be computed |pj|,pO[ . 
This method, which is perhaps less transparent than the one presented in this paper, relates the complexity E(/) to 
the static free energy of r real replicas of the system in the following way: E(r) = (3r 2 ^r and f(r) = a \ ^ Fr ^ ; where 
rF r is the free energy density of the r coupled real replicas in the limit where the coupling goes to zero. To obtain E 
it is therefore necessary to compute F r . For systems with 1-RSB one usually assumes a one step of breaking also for 
the overlap matrix related to F r , with a fixed breaking parameter equal to r and two variational parameters qi and qo 
[^9|,^lJ. In this case the expression of F r as a function of qi, q and r is analogous to the usual static one, except for 
the fact that it has to be minimized with respect to go and gi, but not with respect to r. With these assumptions it 
is possible to obtain an explicit expression of E in terms of /, gi, go and r, where q\ and go are given by their saddle 
point values, while r is fixed by the above equation for f(r). If we now compare this expression for E with the one 
coming from our calculation, we find that they are perfectly consistent: the equations for qi, qo and r are exactly the 
same as the ones for s%, sq and (3lo, and E(r) coincides with jl5|). 
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IV. BEHAVIOR OF THE COMPLEXITY IN THE TEMPERATURE AND FIELD PLANE 



From equations (Q, (|T^) and we obtain the explicit expression of the complexity 

— / »% V n-1 h 2 1, 1, 1 U)Sn 

£(/) = Jwtl«? + yWtl - - log*! + - logZ + -— ^ 

+^ (/ - + ^ log(l - -i)) + ^ (*? - 4) + y " so), (25) 

with t\ = 0(1 — si) and si, so and w given by the saddle point equations (j20|). 

Let us now analyze the behavior of the complexity in the (T, h) plane. For small enough values of temperature 
and field many pure states are present, therefore we expect the complexity to be different from zero in a finite range 
of free energy densities. Actually this is what happens. In Fig. 1 we show £(/) as a function of / at fixed small 
values of T and h: £ is defined between a minimum free energy density fo where it is zero, and a threshold free 
energy density fth where it takes its maximum value. Since for high values of temperature and field only one (finite 
magnetization) paramagnetic state is present ||, we expect that increasing T and h the number of existing states 
progressively decreases until one single state is left. Indeed, if we look at £(/) at a fixed h, but at different values of T, 
we find that the interval of free energies [fo, fth] where £ is defined, becomes smaller and smaller as the temperature is 
increased, and finally it shrinks to a single point, with £ = 0, at a certain critical temperature T c (h). For T > T c (h), 
£(/) is defined only in one point, where £ = 0, and / here coincides with the replica symmetric static free energy, 
that is, the free energy of the single paramagnetic state. In Fig. 2 the line T c (h) is shown in the (T, h) plane p2[ . 
This line separates the region of the plane where only one pure state exists (T > T c (h)) from the region where many 
pure states are present (T < T c (h)) and therefore we identify it as the geometrical transition line. We note that this 
line is monotonically decreasing with the field, so that there is no 're-entrance' p3| . We find that the geometrical 
transition line coincides with the line where the RSB static solution ceases to exist 0. This indicates that, as in 
the case with h = the states with the lowest free energy density are the last to disappear. We note that the 
critical line T c {h) does not coincide in the whole plane with the transition line given by the statics (dotted line in 
Fig. 2) Qj. Indeed the two lines coincide only for h > h* , where the static transition is continuous (qi — qo at 
the transition), while they are strictly different for h < h* (discontinuous transition) ||. This is consistent with the 
different physical scenarios corresponding to the continuous and the discontinuous transitions. In the first case, the 
transition corresponds to different states which merge into a unique ergodic component, so above the transition line 
only one state is present. In the second case, many pure states with finite complexity exist both above and below 
the transition line, and the transition corresponds to the point where the lowest states (/ = /o = /rsb) become the 
relevant ones from a thermodynamical point of view HU • 



V. THRESHOLD ENERGY AND DYNAMICAL BEHAVIOR 



As mentioned in the previous section, at fixed temperature and field, the complexity is defined within a certain 
range [fo, fth] of free energy densities. It is well known that at h = in the low temperature region the states with 
/ = fth ar e the dynamically relevant ones, that is, their energy density coincides with the asymptotic value of the 
dynamical energy ^ n this section we show that this also holds at non zero values of the field. To check this point 
we will perform a dynamical analysis of the model with h ^ in order to compare the dynamical asymptotic energy 
with the energy of the threshold states of the complexity. The equivalence between the threshold and dynamical 
energies can be confirmed numerically for any temperature and field, but for the sake of simplicity we will only quote 
the analytical results at zero temperature and small field. If we take the limit T — > of Eqs. (^0|) and expand it for 
small fields we get the following expression for the energy density of the threshold states, 



'2(p-l) {, , (P-1), 2 \ , nr^y 



ffth = -y y p ' (i + ^y^ h ) +°( h )> T ^°- ( 26 ) 

Let us now compute the asymptotic dynamical energy ]3| (for a more complete dynamical analysis in similar cases 
see [Q,^3|). The relaxational dynamics for the system is given by the Langevin equations, 

SH 

Ht) = -Pj^ - y(t>i{t) + Vi{t), (27) 



G 



where y(t) is a Lagrange multiplier which enforces the spherical constraint, and fi(t) is a Gaussian noise with zero mean 
and variance 2. The dynamics is completely determined by the equations for the two-time correlation and response 

functions, C(t,t') — ± J2i i s i(t)si(t')) and G(t,t') = ^J2i dhtlt') ' and for tne magnetization m(t) = -^ £V ( s »(*))- 
These equations are P,M, 



d t C(t,t') = -y(t)C(t,t') + fi / <M'C p - l {tJ')Gtf,t") 

Jo 

+A»(p-1) / dt"G{t,t")C p - 2 (t,t")C(t",t') + f3hm(t'), (28) 
Jo 

5 t G(t, = -y{t)G(t, t') + S(t - t') + p(p - 1) y <ft"G(i, t")C p - 2 (t, t")G(t", t'), (29) 
m(t) = -y(t)m(t) + nip - 1) / dt"G(t,t")C p - 2 (t,t")m(t")+ [3h, (30) 



where /i = pf3 2 /2. Another equation for y(t) is obtained self-consistently exploiting the spherical constraint, 

y(t) = l-p(3£(t)-(p-l)(3hm(t), (31) 



where £(t) = N 1 (H(t)) is the dynamical energy density. In the limit t — > oo we have the following expression for 
the asymptotic dynamical energy 8^ — lim t _>oo £ (t), 

£oo = - — Jr^ — - — -hrrioo, (32) 
VP P 

where yoc and moo are the corresponding asymptotic values for the Lagrange multiplier and the magnetization. In 
order to compute £oo we then have to find y x and by solving asymptotically the set of dynamical equations 

In ^ it has been shown that for large times there are two regimes for the correlation and response functions. 
The first regime corresponds to time separations r = (t — t') such that r/t — > 0, where the equilibrium fluctuation 
dissipation theorem (FDT) holds, so that Gfbt(t) — — <9 t Cfdt(t)- The second regime, known as aging regime, 
corresponds to r/t ~ O(l). Here time translation invariance is violated and FDT cannot be applied. The response 
function is related to the correlation function by G ag (t,t') — Xddt>C ag (t,t'), where Xd parameterizes the violation of 
FDT. The asymptotic values for the correlation functions in both regimes are, 

lim C F dt(t) = lim C ag (t,t') = q, lim C ag (t, t') = q . (33) 

T^OO t'/t~>l t'/t-tO 

In the limit r — > oo, the equation for Cfdt^) yields a relation between y^ and q which reads, 

y 00 = (l-q)- 1 +Kl-q p - 1 ) ■ (34) 
If we take the limit t'/t — > 1 in the equation for G ag {t, t') we obtain the equation for q, 

^(l-qf-ip-iy^O. (35) 

It is now clear that neither q nor y x depend on the field. Therefore, all the dependence of the energy on the field is 
given by the second term of Eq. (|3^) , while the first term corresponds to the asymptotic energy in zero field ]3j . 

To obtain the asymptotic magnetization m,oo we need to solve the coupled equations for m,oo, qo and Xd coming 
from Eq. (|30|), and from the limits t'/t — > and t'/t — > 1 of the equation for C ag (i, t'), 

= -j/ooTOoo +M [l - <7 P_1 (1 - a;)] +/3h, (36) 

= -yoc9o + - ?)?o 1 + <"9o(l - 9 P_1 ) + l-LXdiq^ 1 - q^ 1 ) + f3hm oc , (37) 

= -yoo + 1 + fi(x - l)q p + m(1 - ) +f3h moo . (38) 
The analytic solution for small /i has a simple form, 

moo = (p-l)(l-q)ph + 0{h 2 ), (39) 
9o = 0(h 2 ), (40) 
Xd= ( P -2)(l- 9) +0(fe2) _ (4i) 
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From here we obtain the dynamical energy, 



(p-2)(l- 9 ) 



(p-l) 2 (l-q), 



P 



0(h 3 ). 



(42) 



In the limit T — > it can be seen from equations (J35|), ( |42[ ) and ( |26| ) that 5^ = E t h- As already said, it is possible to 
check numerically that these two energies are the same in the whole plane (T, h). 



VI. CONCLUSIONS 



In this paper we have computed the complexity E of the states in a very general 1-RSB mean-field model for 
disordered systems, namely, the p-spin spherical model in presence of an external magnetic field. 

We find a region of the (T, h) plane where E is a monotonic increasing function of the free energy density /, defined 
in an interval [/q, /th]- The lower band edge /o, defined by E(/ ) = 0, is always equal to the free energy density of the 
static 1-RSB solution, that is, /o = /rsb- On the other hand, the threshold value f t h gives the free energy density of 
the states asymptotically reached by the non-equilibrium dynamical evolution of the system. 

Our most important result is the one to one mapping we are able to find between the equations and the parameters 
of the complexity at / = fo and the ones of the 1-RSB static approach. This shows that the two calculations are 
essentially the same, even if the ansatz taken for the overlap matrices seem so different in the two cases. 

We would like to add some further comments on what is the main interest of the present work. There are different 
methods by which the complexity of a model can be computed pi|p0| , [l9| . The one used in this paper is the most 
intuitive and transparent, since it simply counts the number of minima of the mean-field free energy function. However, 
from a technical point of view, this method is far from being straightforward and up to now a clear quenched calculation 
of the complexity was lacking. Many efforts have been done in this direction for the Sherrington-Kirkpatrick model 
B, which has a full RSB static solution [p6| p8| , pT| . However, in that case neither the physical meaning of the adopted 
ansatz nor the eventual consistency of the lower band edge with the static free energy density were completely 
satisfying. In the present paper, in the context of a simpler 1-RSB model, we find in a clear way how the computation 
has to be performed and what is the meaning of the ansatz. Therefore, we now believe to have this method under 
control. 

It would be interesting to use this same method for a deeper analysis of the free energy landscape. Indeed, it 
is generally believed that in the dynamics of glassy systems a crucial role is played not only by the minima of the 
Hamiltonian, but also by unstable stationary points ]29| , |T3| ] . In this context, an entropic computation of stationary 
points of any nature is of primary interest. 



ACKNOWLEDGMENTS 



We thank Andrea Crisanti, Leticia Cugliandolo, John Hertz, and especially Jorge Kurchan and Giorgio Parisi for 
helpful discussions. We also thank David Sherrington for useful suggestions and for carefully reading the manuscript. 
The work of A.C. and I.G was supported by EPSRC Grant GR/K97783. The work of J.RG. was supported by EC 
Grant ARG/B7-3011/94/27. 



I " > E-mail: a.cavagnal@physics.ox.ac.uk. 

E-mail: j.garrahanl@physics.ox.ac.uk. On leave from Departamento de Ffsica, Universidad de Buenos Aires, Ciudad 

Universitaria, 1428 Buenos Aires, Argentina. 

E-mail: i.giardinal@physics.ox.ac.uk. 
[1] M. Mezard, G. Parisi and M. A. Virasoro, Spin Glass Theory And Beyond, (World Scientific Publishing, Singapore, 1986) 
[2] A. Crisanti and H.-J. Sommers, Z. Phys. B 87, 341 (1992). 
[3] L.F. Cugliandolo and J. Kurchan, Phys. Rev. Lett. 71, 173 (1993). 
[4] L.F. Cugliandolo and P. Le Doussal, Phys. Rev. E 53, 1525 (1996) 
[5] D. Sherrington and S. Kirkpatrick, Phys. Rev. Lett. 35, 1792 (1975). 
[6] L.F. Cugliandolo and J. Kurchan, J. Phys. A 27, 5749 (1994). 



8 



[7] D.J. Thouless, P.W. Anderson and R.G. Palmer, Philos. Mag. 35, 593 (1977). 

[8] A. Crisanti and H.-J. Sommers, J. Phys. I France 5, 805 (1995). 

[9] J. Kurchan, G. Parisi and M.A. Virasoro, J. Phys. I France 3, 1819 (1993). 
[10] D.J. Gross and M. Mezard, Nucl. Phys. B 240, 431 (1984). 
[11] T.R. Kirkpatrick and D. Thirumalai, Phys. Rev. B 36, 5388 (1987). 
[12] H. Rieger, Phys. Rev B 46 14655, (1992). 

[13] A. Cavagna, I. Giardina and G. Parisi, Phys. Rev B 57 11251, (1998). 

[14] C. Becchi, A. Rouet and R. Stora, Phys. Lett. 52B, 344 (1974); Ann. of Phys. 98, 287 (1976); I.V. Tyutin, Lebedev 

preprint FIAN, 39 (1975). 
[15] G. Parisi and N. Sourlas, Phys. Rev. Lett. 43, 744 (1979); Nucl. Phys. B 206, 321 (1982). 
[16] J. Zinn- Justin, Quantum Field Theory and Critical Phenomena, (Clarendon Press, Oxford, 1989). 
[17] M. Mezard, G. Parisi, N. Sourlas, G. Toulouse and M.A. Virasoro, J. Physique 45, 843 (1984). 
[18] M. Mezard, G. Parisi and M.A. Virasoro, Europhys. Lett. 1, 77 (1986). 
[19] R. Monasson, Phys. Rev. Lett. 75, 2847 (1995). 
[20] S. Franz and G. Parisi, J. Phys. I France 5, 1401 (1995). 
[21] M. Potters and G. Parisi, J. Phys. A 28, 5267 (1995). 

[22] We thank David Sherrington for first suggesting the existence of such a line in the phase h < h* . 
[23] J. Hertz, D. Sherrington and Th.M. Nieuwenhuizen, in preparation. 
[24] H. Sompolinsky and A. Zippelius, Phys. Rev. B 25, 6860 (1982) 
[25] A.J. Bray and M.A. Moore, J. Phys. C 13, L469 (1980). 

[26] C. De Dominicis, M. Gabay, T. Garel and H. Orland, J. Physique 41, 923 (1980). 
[27] S.A. Roberts, J. Phys. C 14, 3015 (1981). 

[28] A.J. Bray and M.A. Moore, J. Phys. C 13, L907 (1980); J. Phys. A 14, L377 (1981). 
[29] J. Kurchan and C. Laluox, J. Phys. A 29, 1929 (1996). 



9 



0.010 



0.008 



0.006 



0.004 



0.002 



0.000 — 
—0.835 




FIG. 1. The complexity £(/) as a function of the free energy density / for T = 0.2, h = 0.2 and p = 3. The minimum value 
S = occurs at /o = /rsb = —0.8322. The maximum value is at f t h = —0.8213, which is the free energy density of the states 
reached by the dynamics. 
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FIG. 2. The full line corresponds to the critical temperature line T c (h) where E reduces to a point and the RSB static 
solution disappears, for p = 3. The dotted line corresponds to the static transition line. The two lines are coincident for fields 
h > h* = 0.408. The static transition is continous for h > h* , and dicontinous for h < h* . The region C is the region of 
coexistence of RS and RSB solutions. 
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